Perfect quantum state transfer with spinor bosons on weighted graphs 
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A duality between the properties of many spinor bosons on a regular lattice and those of a single 
particle on a weighted graph reveals that a quantum particle can traverse an infinite hierarchy 
of networks with perfect probability in polynomial time, even as the number of nodes increases 
exponentially. The one-dimensional 'quantum wire' and the hypercube are special cases in this 
construction, where the number of spin degrees of freedom is equal to one and the number of 
particles, respectively. An implementation of near-perfect quantum state transfer across a weighted 
parallelepiped with ultracold atoms in optical lattices is discussed. 
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The quantum walk (QW) is the quantum mechanical 
extension of the classical random walk, where the quan- 
tum particle (the walker) can be thought of as follow- 
ing many trajectories simultaneously 1]. For quantum 
walks on regular (unweighted) linear graphs, the mean 
displacement of the particle is quadratically faster than 
that of a classical random walk |3| • The time needed to 
propagate from the input to the output vertex for cer- 
tain two-dimensional regular graphs with randomization 
has been shown to be exponentially faster than with any 
classical algorithm , but like the linear graph the prob- 
ability of hitting the output vertex decreases polynomi- 
ally with the graph width. Perfect output probability can 
be achieved on the n-dimensional hypercube 0, 01 based 
on regular and weighted one-dimensional graphs [g with 
an exponential speed-up over a classical random walk, 
though the improvement is only polynomial compared 
with a suitable classical algorithm |^. 

Quantum walks on weighted graphs have been pro- 
posed as an efficient way to transfer quantum states 
(and therefore quantum information) with perfect fidelity 
without requiring external control f^. To maximize 
throughput, it would be most convenient to transport 
many quantum states simultaneously. Unfortunately, 
the quantum particles on which the information is en- 
coded are usually indistinguishable; if the receiver needs 
to know both the information and the particular carrier 
(where the latter might represent a given wire in a quan- 
tum circuit, for example), then the particles need to be 
spatially or temporally separated. For many applications 
it is sufficient to know only that a certain amount of in- 
formation has been transmitted, however. For example, 
Alice might want to send Bob M copies of a given quan- 
tum state (entangled or otherwise), or one qudit each of 
M Bell pairs so that Alice and Bob can share M ebits, 
or — 1 qudits of a GHZat or complete graph state, etc. 

As discussed in detail below, many spinor bosons un- 
dergoing continuous-time quantum walks (CTQWs) can 
be mapped to a single particle walking on an infinite hi- 
erarchy of weighted graphs, including (but not limited 



to) weighted hyperparallelepipeds, hypertetrahedra and 
hyperoctahedra. These graphs generally have vertices of 
variable degree, and share the property with the hyper- 
cubes that they can be traversed with unit probability 
in polynomial time. Unlike the hypercubes, however, 
there is no known classical algorithm that can accom- 
plish the same task. The unique properties of walks on 
these graphs may aid in the development of new schemes 
for quantum communication and computation. 

Consider N bosons with S spin (or pseudospin) degrees 
of freedom, located on a connected graph with V vertices 
of minimal degree one. The {N„ + V - iy./NJ/{V - 1)! 
accessible quantum states for each spin projection a G 

are defined by the 
Bose occupation of the various sites; that is, the natu- 
ral basis states are defined in Fock (occupation number) 
space. Explicitly, for a two-site lattice the HctI-^o- + 1) 
states are 
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where the quantum numbers n.i^ G {0, 1, 2, . . . , A^} de- 
note the number of bosons with spin projection a in the 
site i, |0) denotes the particle vacuum, and X]i<j '^io' ~ 
In this second-quantized notation, the operators and 
flj^ create and annihilate, respectively, a boson with spin 
projection a in site i. The normalization factors reflect 
the Uia-l ways that Uia identical bosons can be arranged 
among themselves. Thus each basis state encodes the 
symmetric permutation of all quantum states with the 
same Hamming weight. The basis states for A^ bosons 
on V sites are obvious generalizations of these, requiring 
at most VS quantum numbers riio-, ^20-, . ■ . , nya-- 

The dynamics of this system are described by the hop- 
ping Hamiltonian: 
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where Ty is amplitude to hop between sites i and j, and 
the angle brackets denote a sum over only nearest neigh- 
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bors of a given vertex connected by an edge. The hopping 
amphtudes therefore define the adjacency matrix of the 
graph, which is assumed to be undirected (r^ = Tji) but 
otherwise unconstrained. The only non-zero matrix ele- 
ments of the Hamiltonian ^ connect an iV-particle Fock 
basis state to another where one boson with spin projec- 
tion a is annihilated at site i and another with the same 
spin is created in an adjacent site j. Consider, for exam- 
ple, the matrix element connecting the states 712, 713) 
and 712 — 1, 713-1- 1), assuming all spin states are equal. 
Because each quantum number labels a different (site) 
subspace, the matrix element decomposes into a tensor 
product of matrix elements for the three sites. If the 
hopping amplitude between these states is r one obtains: 

(tii, 7i2 - 1, 7i3 + l\H\ni,n2,n3) 
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All other elements of the Hamiltonian can be obtained 
by suitably replacing any of the Ui. 

The A'^-particle quantum walk on the original (pri- 
mary) graph is therefore dual to a one-particle walk on a 
secondary graph whose vertices are labeled by the Fock 
states and whose adjacency matrix is defined by the 
Hamiltonian matrix elements above. Since the Hamil- 
tonian is constant in time, the CTQW walk is effected 
by propagating the Fock-space wavefunction according 
to \tp{t)) = U{t)\^{0)) where the time-dependent unitary 
is U{t) = e~*^*. Starting with unit probability on the in- 
put vertex, the CTQW is executed until the probability 
at the output vertex is maximized; this defines the hitting 
time th- The choice of input and output vertices depends 
on the structure of the graph, and will be clarified below. 

Suppose that the primary graph P is the two-vertex 
chain, with the input and output vertices correspond- 
ing to left and right sites, respectively; this can be tra- 
versed with 100% probability with a CTQW in constant 
time th = (7r/2)r-i 0. This value of th is indepen- 
dent of N because each boson executes its own indepen- 
dent quantum walk on the primary graph. The origi- 
nal eigenvalues A — {^t} are transformed into the A^- 
muhiples A = t{-N, ~N + 2, . . . , N - 2, N} which are 
always commensurate jg] . The resulting structure of the 
evolution operator U (t) on the secondary graph ensures 
that perfect state transfer occurs between each pair of 
symmetry-related vertices a,tt — th, i.e. between any two 
vertices whose quantum numbers are related by inversion 
\nia, n2a, ■ ■ ■ , nva) ^ \nv<T, ■ ■ ■ n2a, nia). 



Choosing N = S so that there are N distinguishable 
bosons, then riia € {0, 1} for each value of a and the 
matrix elements ^ are all equal to t. As expected, 
the system maps to the A-dimensional r-length hyper- 
cube P'^^ ^ which has the same value of th- When 
5 = 1 so that all the bosons are spinless and there- 
fore indistinguishable, the secondary graph Gn is the 
weighed line whose A -I- 1 vertices are labeled by the 
Fock states | A — 71, ti) where n G {0, 1, . . . , A^} and the 
edge weights are given by Hn,n+i — r-y/ (n + 1)(A — n). 
These are precisely the weights required to transport 
a quantum state across a quantum wire, also in time 
th — (''■/2)t~-^ 0] which has been shown to be optimal 0|. 
In fact, as has been noted these weights result from 
projecting the vertices of an A^-dimensional unit-length 
hypercube onto a line. The boson construction indicates 
that the projection is driven by indistinguishability: all 
states of the hypercube with the same Hamming weight 
are mapped to a common point on the line, but the bit 
positions in the original strings are not specified. 

Suppose that the A bosons each carry the qudit state 
ItA) = cil- (^) + . . . + C5|^), where Ef=i IcP = 1. 
Because the Hamiltonian Q preserves the total spin, 
the A'^-particle walk decomposes into a superposition of 
walks on S graphs of the type Gn (quantum wires), 
S{S — 1) graphs of type Gn~i P (weighted rectan- 
gles), S{S - 1){S - 2)/2 graphs of type Gn-2 O 
S{S — 1) graphs of type Gn-2 d) G2, etc. These sec- 
ondary graphs are hyperparallelepipeds of maximum di- 
mension S, where along each direction the edge weights 
correspond to those of the quantum wire. Perfect quan- 
tum state transfer is effected between any two vertices 
connected by a line through the geometric center of the 
hyperparallelepiped. 

Consider now A' bosons randomly distributed on one 
of the weighted secondary graphs. At t = th, the quan- 
tum state labeled by the Bose occupations for each site 
is transformed into another state with the same popula- 
tions in the symmetry-related sites. This final quantum 
state corresponds to a site that is itself the symmetry 
counterpart to the initial site on a tertiary graph. In this 
manner one can generate an infinite hierarchy of weighted 
graphs that all share the properties of the primary graph. 
The largest edge weight is proportional to A^^ • • • A^2A^i 
where Nm is the number of bosons at recursion level m; 
rescaling this weight back to r yields a physical hitting 
time th oc YViLi that increases polynomially even as 
the number of vertices grows exponentially. 

The smallest tertiary graphs built from quantum wires 
and hypercubes are shown in Fig. for S" = 1. With 
A^' indistinguishable bosons on Gn, the resulting graphs 
are A^-dimensional weighted hypertetrahedra with A^' -I- 1 
vertices along each edge, but where not all vertices share 
an edge with their nearest neighbors. For the two- 
dimensional graphs defined by A' bosons on G2 , shown in 
Figs.^a-c), the edges between V vertices along a given 
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FIG. 1: Examples of small tertiary graphs built from the 
quantum wire and hypercube. Cases (a) through (c) corre- 
spond to A'^ = 1 through 3 on G2; cases (d) and (e) depict 
N = 2 on G3 (tetrahedron of length 2) and on P®P (octahe- 
dron of base length 2), respectively. All weights are in units 
of T. A CTQW running for t = {tv/2)t~^ will yield perfect 
quantum state transfer between vertices 'in' and 'out'. 

constant direction have weights defined by ^ — 1 bosons 
on P, multipUed by v^; perfect quantum transfer occurs 
between any two vertices related by a reflection about 
the graph center, so that states on central vertices re- 
main fixed at intervals of t^. N' bosons on the square 
□ (two-dimensional unweighted hypercube) map to a sin- 
gle boson on (generally truncated) octahedra with iV' + 1 
vertices along each edge of the square base, as shown in 
Fig. n^e) ; quantum states are transferred to between ver- 
tices through the body center and parallel to the base, so 
that states on the two apices remain fixed at th intervals. 

The edge weights on high-order graphs in this con- 
struction generally increase both with N' and the num- 
ber of vertices, as is evident in Figs, na)-(c). Consider 
for example the graphs defined by a hierarchy of spinless 
bosons {Nm on [. . . Ni on (Nq on 2'*'')]}, where d is the 
secondary hypercube dimension and m is the number of 
iterations in the hierarchy used to construct the graph. 
Because of the symmetry of the underlying graphs at each 
level, edges with the largest weights are contiguous. This 
implies an efficient classical algorithm to traverse these 
graphs with high probability and relatively few steps: one 
simply chooses the 'heaviest' edges at each step. 

Graphs avoiding this issue include constructions such 
as {Nm on ... [Ni on {Nq on 2)f^' . . .}»^", with iVi 
distinguishable bosons on hypertetrahedra forming the 
underlying graph for N2 distinguishable bosons, etc. An 
example is shown in Fig. El When these graphs get very 
large, the weights tend to become homogeneously dis- 
tributed; because the vertex degrees also vary across the 
network, there are no identifiable criteria that can be 




FIG. 2: Two representations of the same four-dimensional 
fourth-order graph (2 on G2)®^ = [2 on (2 on 2)]®^ Each 
triangular node in (a) corresponds to the six-vertex weighted 
graph shown in Fig.^b); each node within a triangle connects 
to its counterpart in an adjacent triangle by an edge with the 
same weight. Solid and striped lines correspond to edges with 
weight y/2 and \/4, respectively. 

used to classically determine the optimal edge at each 
step. Together with an exponential increase in the num- 
ber of vertices, this implies an exponential speed-up of 
the quantum walk over any classical traversal algorithm. 

The two-vertex graph need not be the only primary 
graph that yields perfect quantum state transfer, how- 
ever, if the weights are allowed to be complex. The re- 
sulting hierarchy of higher-order graphs are the quantum 
analog of signed classical networks. In this case the hop- 
ping Hamiltonian becomes H = ^io-^icr +H-C.), 
where r is now complex and 'H.c' stands for 'Hermitian 
conjugate.' One of these is the triangle A with edge 
weights ir/V^ and eigenvalues A — {±t, 0}. Labefing 
the vertices clockwise around A, the state is transfered 
to each successive neighbour in a counterclockwise fash- 
ion after each time interval t = {2tt/3)t~^ . The sec- 
ondary graphs generated with N bosons on A corre- 
spond to weighted triangles and tetrahedra. These are 
much like Figs.^a-c), in that the weights are defined by 
the Gv graph with V + 1 vertices along a given direc- 
tion (multiplied by j/a/S), except that all nearest neigh- 
bors now share an edge. Quantum states are successively 
transferred between vertices related by a counterclock- 
wise 27r/3 rotation about the center, which is a symmetry 
operation for the graph. Another example is the square 
□ or bowtie ixi with successive weights 1 , e'^°' , e^^ , e^'^ in 
units of T, where the first weight can be chosen to be real 
without loss of generality; P (g) P is a special case with 
a — P = "/ = 0. Whenever 7 = a + f3, the eigenvalues 
are A — {±t} and perfect state transfer occurs between 
vertices at opposite corners. 

That the eigenvalues of the hopping Hamiltonian on 
the primary graph are commensurate and symmetric 
about zero is a necessary but insufficient condition for 
perfect state transfer, however. For example, the com- 
plete four-vertex graph can have weights a on outside 
edges and weights b on the inside (diagonal) edges, while 
still preserving 7r/4-rotation symmetry. For example, 
choosing a = i/y/2 and b = i yields A = {±2r, 0} but 
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after th = 7t/2t the state is in an equal superposition 
of two vertices. While it is conceivable that a judicious 
choice of weights could guarantee perfect state transfer 
between two nodes of an arbitrary graph, the prescription 
for accomplishing this has not been determined. 

Ultracold atomic spinor bosons confined in optical lat- 
tices 0, 0, provide a physical system that is particu- 
larly amenable to implementing perfect quantum state 
transfer on some of the weighted networks discussed 
above, especially those with regular geometries. The 
atomic interactions can be made small through the use of 
Feshbach resonances ^^'^ ^^e tunneling amplitudes 
and phases from site to site can be controlled to some 
extent |l2l | . Particularly simple graphs to physically con- 
struct are the secondary weighted parallelepipeds G®'' 
where 1 < d < 3 and the weighted tetrahedra based 
on A with complex weights. The tunneling matrix ele- 
ments along any given direction now depend on position, 
Ti_i^i = Ty/i{V - i + 1) with 1 < i < V. With the 
substitution i = [(V" + l)/2] — j where j is the vertex 
index relative to the graph center, one finds tJ^ij oc 
[{V + l)/2]^ — j^. For standard optical lattice potentials 
of the form V{x) = sEncos^^kx), where x = {x,y,z}, 
k = 2tt/X is the wavevector (A is the laser wavelength), 
and s is the lattice depth in units of the atomic re- 
coil energy En = ?i^fc^/2m (m is the atomic mass), the 
tight-binding approximation applicable for s > 5 gives 
= (16/7r)s'^/^e~'*^ in recoil energies The re- 

quired quadratic variation ofTj_i j can be approximately 
obtained by choosing a weak position-dependent lattice 
depth s — sq + ap in each direction, which in principle 
can be accomplished by appropriately focusing the lat- 
tice laser beam and adding end caps to provide hard-wall 
boundary conditions. 

The Schrodinger equation for bosons in a one- 
dimensional (ID) optical lattice with V — 101 and 
So = 5 was propagated in time numerically using a 
finite-element discrete variable approach Choosing 
a = 10"'^ and S* = 1, an initial quantum state centered 
s-t jo = {10, 20, 30, 40} and spread over several sites was 
found to transfer to the sites located symmetrically oppo- 
site the lattice center with probabilities exceeding 99%. 
The hitting times were within 10% of the theoretical 
value th = {V + l)(7r3/Vl6)s-3/4e2v^(^/^^) ^ 40 ms, 

assuming ^'^Rb atoms in a lattice with A = 800 nm. Sur- 
prisingly, the initial wavepacket remains well-localized 
during the propagation. Weak damping of the oscilla- 
tions between initial and final states is due to increased 
dispersion driven by the imperfect phase profile (which is 
more pronounced near the lattice edge), and was found 
to have a characteristic time between A2th (jo — 40) and 

iietft (jo = 10). 

Precise control of the lattice potential or atomic inter- 
actions is not required to effect the state transfer. For 
example, choosing a — 7 ■ 10^^ when V = 101 gives a 
variation of the hopping amplitudes that is more appro- 



priate for a lattice with V = 41. State transfer with 
probability exceeding 95% was found for any jo < 20 
within approximately 10% of th — 380t^^. For jo > 20, 
the wavepacket remained at the lattice periphery. Vir- 
tually identical behavior was obtained using a Gaussian, 
rather than quadratic, profile for the lattice height. Weak 
repulsive interactions between bosons were incorporated 
in a mean-field sense by solving the appropriate quasi- 
ID nonlinear Schrodinger equation. The strength of the 
additional term g\ip\'^ is parametrized by the effective ID 
interaction strength 17, which in the current dimensionless 
units is g = asX/'!T'^i± 0| where Og is the s-wave scatter- 
ing length and £± is the transverse oscillator length. For 
parameters relevant to current experiments (1(1 |. g « 0.2. 
For values of 5 < 1, t/i is found numerically to increase 
only slightly (by 3%), and the probability of hitting the 
output vertex decreases to 95%. Because the th are iden- 
tical for all sites, the behavior for this range of interaction 
strengths is only weakly dependent on the initial density 
distribution. For larger values of g the output proba- 
bility drops precipitously, but the mean-field approxima- 
tion becomes increasingly suspect. The numerical results 
suggest that near-perfect quantum state transfer can be 
achieved using ultracold bosons in weighted optical lat- 
tices under realistic experimental conditions. 
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